Electron-phonon coupling and longitudinal mechanical-mode cooling in a metallic 

nanowire 
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We investigate electron-phonon coupling in a narrow suspended metallic wire, in which the phonon 
modes are restricted to one dimension but the electrons behave three-dimensionally. Explicit the- 
oretical results related to the known bulk properties are derived. We find out that longitudinal 
vibration modes can be cooled by electronic tunnel refrigeration far below the bath temperature 
provided the mechanical quality factors of the modes are sufficiently high. The obtained results 
apply to feasible experimental configurations. 
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Electron-phonon coupling in metals, albeit extensively 
studied over several decades [l( , is of utmost interest and 
importance in view of present day developments in nano- 
electromechanics J3j3| and in electronic cooling and sens- 
ing on nanoscale [J, [5| . A number of questions arise when 
the dimensionality of the phonons is reduced from the 
conventional bulk three-dimensional case [f| 0, H] . Re- 
cent experimental observations of metallic wires on thin 
dielectric membranes support the fact that reduction of 
phonon dimensionality leads to weaker temperature de- 
pendence of the heat flux between electrons and phonons 
[9(. Very little is known about truly one-dimensional 
wires, where transverse dimensions are far smaller than 
the thermal wavelength of the phonons, although this 
regime is readily available experimentally at sub-kelvin 
temperatures in wires whose diameter is of the order of 
100 nm or less. Recently though, substantial overheating 
was conjectured to be the origin of excess low-frequency 
charge noise in a suspended single-electron transistor in 
this particular one-dimensional geometry [10| . In this 
Letter we derive an explicit result for electron-phonon 
heat flux in a metallic wire in which electrons behave 
three-dimensionally but phonons are confined to one di- 
mension, and relate this result to the standard bulk result 
for the corresponding metal. We present a scenario of 
tunnel coupling the metal electrons in a wire to a super- 
conductor on bulk, whereby cooling of wire electrons can 
be realized. We demonstrate that the few available me- 
chanical modes, i.e., discrete longitudinal phonons, can 
be cooled significantly by their coupling to the cold elec- 
trons in the wire. This occurs provided the mechanical 
modes are not too strongly coupled to the thermal bath, 
meaning that the mechanical Q value of the mode is suf- 
ficiently high. Recently, indirect experimental evidence 
of electronic coo ling of phonons in a bulk system was put 
forward in Ref. [111 ]. 

To obtain results for the electron-phonon heat flux in 
a one-dimensional metallic wire (see Fig. [T]for the geom- 
etry and thermal model), we follow the standard proce- 
dure from the existing literature normally applied to the 
case of either bulk three-dimensional phonons [l], [Tl] , or 
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FIG. 1: The system under study. In (a) we show the sus- 
pended wire whose transverse dimensions are supposed to be 
smaller than the thermal wavelength of the phonons, Athermai- 
In this particular example the metal wire is connected to the 
bulk superconducting reservoirs via tunnel barriers to form 
a tunnel junction refrigerator. In (b) we show the relevant 
thermal model of the system. 



to the case where phonons are restricted to a semi-infinite 
bulk [f|. The net heat flux from electrons into a discrete 
phonon mode [i at wave vector q is given by 

Q e ^(q) = 2 £ ^[r* (k -> k - q ) - r^(k -> k + q )], 

k 

w 

where phonon emission (e) and absorption (a) rates by 
the electrons with wave vector k are obtained via the 
golden rule as 

r.(^ k - q ) = f| v ,|>(^-) + i] 

«/M - /(Ek-,)]J(E k - £ k _„ - fiu„), (2) 



and 
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x/(£ k )[l - /(£ k+q )]rJ(£ k - £ k+q + (3) 
Here g M and n{j£*r) = [exp(^-) - 1]" 



are the 
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electron-phonon coupling constant and the Bose distri 
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bution, respectively, of the phonon mode p at angu- 
lar frequency uj^ and at temperature T M , and f{E) = 
[exp( fc ^ T ) + is the Fermi distribution of the elec- 
trons at temperature T e . 

We next evaluate g^q using the standard results of 
deformation potential of the collective lattice vibrations 
pj. Let w M (r) be the displacement vector for mode 
\x normalized in the volume of the wire, V, such that 
/ d 3 r w M (r)-w*, (r) = 5^'. Then can be obtained 
from the divergence of w M (r) 



9w = l E F)J^J d3r *k-q(')*kW V-w M (r), (4) 

where ^(r) are the electronic wave functions. Here Ep 
is the Fermi energy of the electrons and p is the mass 
density of the wire. 

The thermal wavelength of phonons, Athcrmai = 
at temperature T and mode velocity C£ is typically of 
order 1 pva at T = 100 mK. In a wire whose length 
L > Athcrmai and with transverse dimensions <C Athermab 
only modes with q = (0, 0, q) directed along the wire (z- 
axis) appear relevant, since the ones with perpendicular 
q are too high in energy. There are basically four types 
of vibrations: longitudinal, flexural (two, with x and y 
polarizations), torsional and shear modes [ll] . The last 
one has a gap and is therefore not excited at low tempera- 
tures. Of the remaining ones the torsional modes have no 
divergence, and essentially only the longitudinal modes 
couple to electrons in the long wavelength limit. Exper- 
imentally this seems to be the case in carbon nanotubes 

We consider longitudinal modes with specific boundary 
conditions: the wire (or the three-dimensional body) is 
assumed to be clamped at the ends. As we will detail 
below, this corresponds to a feasible realization. Let the 
wire extend from z = to z = L. Then, the normalized 
longitudinal eigenmodes of the beam are given by 



w p» 



: sin(/i7rz/L)z, fj, = 1, 2, 3, 



(5) 



They are characterized by the linear dispersion relation 
u>n = cipn/L, where eg = y/E/p is the longitudinal 
sound velocity (E is Young's modulus). Assuming zero 
electronic boundary conditions along with equal elec- 
tronic and phononic volumes we obtain again in the long 
wavelength limit 
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where A4"j 



The momentum q transferred 



1 hE% 
9 pVc e ' 

between the electron and the vibrational modes of a 
clamped beam takes discrete values q^ only and is by 
convention positive. 

We perform next the integration over electron ener- 
gies in Eqs. and and insert the results in Eq. (TT|) 



obtaining 

• 2irMjcjmN(E F ) 2 Miq^ 



hkp 



(7) 

Here, m is the electron mass, kp the Fermi wave vector, 
and N(Ep) the electronic density of states at the Fermi 
energy. Three-dimensional distribution of electrons was 
assumed here, since we discuss only the case of ordinary 
metals, where kp 1 <C 1 nm. Using the definition of M. 2 

2 

£ of the free 
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above, and N(Ep) = and E F = 

electron gas, the prefactor in Eq. ([7]) can also be written 



in the form 



2irMfcjmN(E F ) 



h k F C£ 



The total heat 
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flux between electrons and phonons can then be obtained 
as a sum over all modes: 



(8) 



We obtain the continuum result for a long L 3> A t hermai 
ID wire by assuming a uniform density of modes with all 
of them at the same temperature T M = T p . We then 
replace the sum by an integral, ^2 q —f ^ J °° dq. After a 
straightforward integration we obtain 



Here, Sid is given by 
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It is instructive to compare this result to the celebrated 
result for longitudinal phonons in three dimensions (see, 
e -g-j [13 an d references therein) 

Q e ^p = SV(T e 5 - T p 5 )- (11) 
Here, the material specific prefactor S is given by: 



C(5) k F kg 
3tt 3 H 3 cjp 



(12) 



We conclude that Eid is related to the known E of the 
bulk by 



_n((3) hc e 2 

^id — ^tttvIi — I ^- 
6 C(5) k B 



(13) 



Note that Eq. © with the relation (fT3")) between Ei£> 
and E are quite general, and do not depend on the 
choice of free electron gas parameters that lead to 
Eqs. CE2) and flTUJ). Equation © with the help of (JT3J) 
and the experimentally determined E can then be used 
to assess electron-phonon coupling in one-dimensional 
wires. Equation (p~2|) predicts the behaviour of real 
metals rather well: the overall magnitude of E from 
(JT2J with parameters of usual metals is of order E ~ 
10 8 WK~ 5 m~ 3 , whereas measured values are typically 
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around 10 9 WK~ 5 m -3 . The deviation may be partly as- 
cribed to the complicated structure of the Fermi surface 
in real metals Q. 

Equations (fTTj) and @ predict correctly the cross- 
over between three-dimensional and one-dimensional be- 
haviour. To see this, let us look at the linearized heat 
conductance for a small temperature difference AT = 
T e —T p between electrons and phonons, such that Q e ^ p ~ 
G ep AT. From Eq. (JTTJ, we obtain = 5EVT 4 , where 
we denote by T the (almost) common temperature of 
the two subsystems. Similarly from Eq. © we obtain 
= 3Si D iT 2 . Now let us consider a wire whose 
square cross-section iswxro. The cross-over between 3D 
and ID behaviour is expected to occur when the first lon- 
gitudinal modes get occupied thermally within the cross- 
section, i.e., when hci/w ~ fcgT. Making use of the 
relation (|13p. and V = Lw 2 , we then see that with the 
above condition the expressions of G^ and G^ become 
identical in form, apart from numerical prefactors. 

Next we demonstrate that variation of electron temper- 
ature in the wire leads to variation of the temperature 
of its vibrational modes. In particular, electron medi- 
ated mechanical mode cooling becomes possible. If we 
assume a highly underdamped mechanical mode whose 
quality factor 3> 1, we can obtain the heat flux from 
the thermal bath into the mode fj, in a classical picture 
as 




;bath — >p — TT (Tbath t^)' 



(14) 



FIG. 2: Minimum electron temperature in a metal wire as a 
function of bath temperature, cooled by SINIS tunnelling for 
various values of Q. The solid blue lines are from the one- 
dimensional model, and the red dashed lines from the three- 
dimensional model. The parameters we used correspond to Al 
as a superconductor and Cu as the metal wire of L = 1 pm 
length; its width and thickness are assumed to be 30 nm both, 
£ = 2 • 10 9 WKT 5 m- 3 , E = 130 GPa and p = 8920 kgm" :i . 
The tunnel resistances of the two NIS junctions are Rt = 10 
kf2 both, and we assume that the non-ideality parameter of 
the junctions has a value 7 = 1 • 10 -4 (curves lying higher) or 
7=1. 10~ 5 (lower). 



This result can be inferred as a solution of the Fokker- 
Planck equation of Brownian motion in the harmonic po- 
tential or by direct solution of Langevin equation [15l |. 
We have assumed that the mode temperature is given 
by the equipartition principle via ksT^ — k(x 2 } for the 
position x of the Brownian particle with spring constant 
k. Equation (fT4|) is the high temperature limit of the 
quantum expression of heat flux 
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Wa* bath 
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(15) 



which is identical in form with Eq. ([7]). We have again 
identified lu^ = cgq^. One then finds a steady-state tem- 
perature of the mode p, by solving the balance equation 
(see Fig. P, 



Qbath->/i + Qe-^ft — 0. 



(16) 



if H£i. 
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There are some interesting limits: 

° Qn 187T p 

electrons cool efficiently and the mode temperature fol- 
lows T e , whereas in the opposite limit the mode temper- 
ature stays at Tbath- Eliminating &f in favour of exper- 
imentally determined E, we find that the temperature 
of the mechanical mode follows that of the electrons if 
With parameters of ordinary metals 



^> 12 g( 5 ) 



^ Weft, 

this leads to the condition ^> 100. 



We conclude the formal part by obtaining a useful re- 
lation yielding the heat flux between electrons and the 
bath with the help of their respective temperatures, us- 
ing Eqs. 0, ©, (HSJ) and ([15]). and assuming that all 
the relevant modes have the same quality factor Q: 



->bath 



y k B ' 



12C (5) k 



R4 r 3 



x > ,9,>(-r-^r) 



' ksT e 



„ ^ )]• 



k B T h 



ath 



(17) 



An expression of type © can be obtained in the con- 
tinuum limit again, but here the factor Eid must be re- 
placed by [1 + ^^Qq-iy^. 

We next apply the results above to determine the per- 
formance and mechanical mode cooling [l6|, [I?], [H, [H| 
in a suspended electron refrigerator. Note that overheat- 
ing of a suspended wire, or a single-electron transistor 
[1 01 ] , can be analyzed similarly as our example of cooling 
below: heat currents and temperature drops are simply 
inverted. In a hybrid tunnel junction configuration (SI- 
NIS) , with a metal island (N) and superconducting leads 
(S), the electron system in N can be cooled far below 
the bath temperature by applying a bias eV of the order 
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FIG. 3: Population of the lowest vibrational mode in a wire 
against bath temperature according to the one-dimensional 
model. The solid (blue) lines show the population when the 
electron cooling is applied, whereas the dashed (red) lines are 
the corresponding non-cooled populations. The parameters 
are the same as in Fig. [2] The results for the two values of 7 
are almost indistinguishable, and we have plotted only those 
for 7 = 1. 10~ 4 . 



of the superconducting gap A over each tunnel junction 
(I) between S and N. This SINIS refrigeration technique 
based on energy filtering of the tunnelling electrons due 
to the gap in the superconductor has been applied ex- 
tensively over the past decade, for a review see Ref. 
but not yet in suspended wires to the best of our knowl- 
edge. Here we propose its use in connection with the 
one-dimensional phonon system. It is possible to cool 
not only the electrons in the wire but also the vibra- 
tional modes in it by coupling them to the cold elec- 
trons. Figure [2] shows numerically calculated results for 
the minimum electron temperature reached as a func- 
tion of the bath temperature: at the optimum bias volt- 
age of the junctions heat is removed from the wire at 
a rate Q ~ A 2 /(e 2 R T )(T e /T c f /2 . In steady state this 
heat flux is balanced by the heat flux from the phonon 



modes. We assume that all the relevant modes have the 
same quality factor Q n = Q. The collection of results 
in Fig. [2] shows that if Q is large, strong suppression 
of electron temperature can be achieved. The satura- 
tion of the temperature with low Tb a th is caused by the 
ohmic heating in the refrigerating junctions with leak- 
age parameter 7, which has been chosen to correspond 
to typical experimental conditions: 7 equals the low tem- 
perature zero bias conductance of a junction normalized 
by the value of conductance at large voltages, and it can 
be conveniently included in the (normalized) density of 
quasiparticle states of the superconductor at energy E 
as n s (E) = |Rc( , E+nA )| 41. The cooling effect 

of the suspended structure differs from that of the result 
of the three-dimensional model; specifically the results 
of the one-dimensional model, valid when w < A t hermab 
do not depend on the transverse dimensions of the wire, 
whereas the results of the three-dimensional model are 
determined by these dimensions as well via the depen- 
dence on volume in Eq. (fTT|) . Also the vibrational modes 
involved are cooled: this is demonstrated in Fig. [3j where 
we plot the population of the lowest mode, n = 1, un- 
der the same conditions as in Fig. [2l The corresponding 
mode occupations in the absence of electron cooling are 
shown for reference. The magnitude of the mode cool- 
ing is determined by the interplay of the cooling power, 
electron-phonon coupling, and the coupling to the bath, 
determined by Q. From our example it seems obvious 
that electron-mediated cooling of the vibrational modes 
into the quantum limit is a feasible option, manifested 
by the very low mode populations, in particular when Q 
is large. 

In summary, we derived the basic relations governing 
electron-phonon heat transport in narrow metal wires, 
where the electron distribution is three-dimensional and 
the phonon distribution is confined to one dimension. In 
this realistic scenario describing suspended wires made 
of ordinary metals, we find that the heat currents differ 
drastically from those in bulk systems. In particular, we 
demonstrated that the vibrational modes of the wire can 
be cooled significantly by electron refrigeration, provided 
the mechanical Q's of the modes are sufficiently high. 
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